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ABSTRACT 


they  give  are  close  enough. 


The  waveform  Bounding  aporoacn  to  fast  timing 
analysis  of  *OS  VLSI  circuits  is  discussed.  The 
idea  is  to  compute  rigorous  closed-form  expressions 
giving  uppe-  ana  lower  pounds  for  transient  voltage 
wavefoms,  ratner  than  exact  values.  Th;  goal  is 
to  enaole  rapia  computation  without  sacrificing 
user  confidence  in  tne  results. 

I.  Background  and  Objectives 

Existing  aporoacnes  to  timing  analysis  and 
simulation  of  digital  integrated  circuits  fall, 
roughly  sneaking,  into  th'-ee  classes: 

i)  Methods  sucr  as  SPICE2  [1]  ana  ASTAP  [2], 
based  on  essentially  exact  numerical  solution  of 
the  network's  differential  eauations,  are  accurate 
ana  reliaoie.  But  even  with  tne  increase  i r,  soeea 
afforded  cv  the  waveform  relaxation  metnod  [1], 
exact  numerical  solution  is  too  slow  for  tne  needs 
of  tne  VLSI  era. 

ii)  Special :aed  MOS  timing  simulators  like 
MQTIS-C  [A]  ana  SPLICE  [5]  rely  on  table  look. p  o' 
device  cnaracten sties  for  speed,  and  save  addition¬ 
al  time  Dy  terminating  a  Newton-Rapnson  or  similar 
iteration  before  convergence  is  reacned.  SPLICE 
is  in  addition  a  mixed-mode  circuit,  timing  and 
logic  simulator  and  uses  a  selective  trace  algo¬ 
rithm  to  exploit  latency.  In  both  these  programs 
the  termination  of  an  iterative  step  prior  to  con¬ 
vergence  saves  time  at  tne  cost  of  accuracy  and,  in 
some  instances,  of  numerical  stability  [6],  The 
improvement  in  speed  over  SPICE2  is  typically  one 

to  two  orders  of  magnitude  for  SPLICE  [5]  and  about 
two  oroers  of  magnitude  for  MOTIS-C  [7], 

iii)  More  recently,  some  researchers  are  ex¬ 
ploring  an  alternate  aporoacn  to  timing  analysis 
and  simulation  based  on  a  radically  simplified 
electrical  description  of  the  network.  RSIM  [8], 
CRYSTAL  [9],  and  TV  [10,11]  fall  at  the  far  end  of 
the  speed-accuracy  tradeoff  curve  from  SPICE2.  A 
MOSFET  is  typically  represented  in  these  programs  by 
an  extremely  simplified  model:  a  linear  resistor  in 
series  with  a  switch.  And  a  oolysilicon  or  dif¬ 
fusion  line  is  represented  by  a  lumped  capacitance 
in  RSIU,  or  By  a  delay  in  TV  obtained  by  simply 
averaging  the  upper  and  lower  delay  bounds  obtained 
by  Rubinstein,  Penfield,  and  Horowitz  [12].  These 
programs  are  potentially  very  fast  and  have  a  num¬ 
ber  of  attractive  user-oriented  features.  The 
drawback,  of  course,  is  that  there  are  no  absolute 
known  limits  to  the  error  in  their  total  delay 
estimates.  The  user  can  never  be  sure  the  answers 


The  objective  of  the  waveform  bounding 
approach  to  timing  analvsis  and  simulation  is  to 
combine  tne  computational  speed  that  results  from 
avoiding  the  numerical  solution  of  differential 
eauations  with  the  user  confidence  in  tne  result 
tnat  comes  from  rigorous  uncertainty  pounds.  Our 
attack  on  tne  timing  analysis  proolem  is  based  cr, 
a  careful  fundamental  study  of  the  differential 
equations  describing  the  dynamics  of  gates,  pass 
transistors,  interconnect,  ano  the  standard  digital 
circuits  constructed  from  them. 

In  additiw  to  the  MIT  croup  working  on  this 
project,  Marx  Horowitz  [’2,13]  is  currently  com- 
plefng  a  dissertation  on  MOS  timinc  analysis  at 
Stanford. 

II.  Response  Bounds  *cr  Interconnect 
2.11  Linear  Interconnect  Mode's 

This  section  summarizes  the  results  octa-nec 
in  [12].  In  this  work  an  MOS  signal  distribution 
network  as  shown  in  Fin.  1  is  modelled  as  a 
branched  linear  RC  line,  i.e.,  an  RC  tree,  as  m 
Fig.  2. 
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Figure  1.  Typical  MCS  signal-distribution  network. 

The  invert er  is  shown  driving  three  gates. 


Figure  2.  The  linear  RC  tree  shown  above  is  a  model 
for  the  network  of  Fig.  1.  The  voltage 
source  is  a  unit  step  at  time  t  *  0. 

For  any  two  nodes  in  the  network,  R^i  is  defined  as 


2 


thesun  of  the  resistances  along  the  route  consisting  of 
the  intersection  of  the  path  from  tha  input  tonode  t  with 
the  oat"  from  the  incut  to  node  m,  as  illustrated  in  Fig. 
3.  *he  three  time 

weet  > 


Figure  3.  I' lustration  o»  resistance  terms.  ror  this 
network,  «k  1  =  S']  *  R^,  Rjdj  =  R^  ♦  R£  »  R3, 
ano  R,-  ,■  =  Ri  +  r  R5 . 

constants  used  to  derive  response  Bounds  are 
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wnere  tne  surnations  are  taken  over  all  nodes  of  the 
network. _  The  derivation  in  [12]  shows  that  vs(t) 
v<( t)  v<it),  for  all  t  »  0,  where  v,-(t!  is  the  actual 
zero  state  step  response  at  any  terminal  node  i,  and 
tne  oounos  v.,'t!  anc  v;;t)  are  given  by 
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as  illustrate: 
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rigure  4.  Form  of  the  oounos  with  the  distance  from  the 
exact  solution  exaggerated  for  clarity. 

The  time  reouired  to  compute  these  bounds  grows  only 
linearly  with  the  numoer  of  elements  in  the  network. 
Recent  applications  of  this  result  include  [10,11,14,15 
16].  The  ultimate  goal  of  this  portion  of  the  project 
is  to  derive  a  hierareny  of  such  bounds,  permitting  the 
user  to  trade  off  accuracy  for  computation  time. 

2.2)  Nonlinearities  Affecting  Interconnect 

The  linear  circuit  model  in  Fig.  2  fails  to 
incorporate  three  types  of  nonlinearities  present  in 
Fig.  1  or  related  circuits:  the  nonlinear  output 


resistance  of  the  inverter,  the  nonlinear  gate-to- 
cnannei  capacitance  0 (  tne  ‘iOSFE"  loads,  and  the  non¬ 
linear  capacitance  from  any  diffusion  line  to  substrate. 
This  section  describes  recent  work  [17-20]  that  allows 
the  Oounos  for  linear  networks  [12]  '0  oe  applied  to 
RC  lines  incorporating  suen  nonlineanties.  (Furtner 
research  is  needed  for  branched  lines,  i.e.  RC  trees.; 

Using  the  notation  and  sign  conventions  illustrate" 
in  Fig.  5,  the 


Figure  5.  Two-capacitor  example  of  a  nonlinear,  non- 
uniform  RC  1 ine. 


state  eauations  for  any  nonuniforr,  nonlinear  lumped 
RC  line  with  N  capacitors  can  be  written  in  the  forr 


j  j 


1  <  j  1  N,  -'6! 

where  gg q  =  0,  v  j  =  e,  and  the  capacitor  constitutive 
relations  Q;.=  h  - ( v s )  are  continuously  differentiate 
with  C,  ( v  j)  —  h*J(v^)>  0  everywhere.  Re  assume  the 
resistBr  curvesJare  continuously  dif<e'ertiaPle,  strictly 
increasing,  ana  pass  through  the  origin. 


Lemma  1  [19] 


Consider  any  nonlinear,  nonunifc'*"  RC  line,  at  any 
instar.t  during  an  "up"  transition  (i.e.  e  >_  0,  e  ^  3) 
from  eauil ibrium, 


v  j  ( t )  _>  0 ,  v  j  ( t )  <_  e ,  v  ■  ( t )  v  _  „  ;  t ) 

and  v- (t)  >  0,  1  <  j  <  N. 

Lemma  1  is  proved  in  [19],  Using  it,  we  give  1 
proof  in  [20]  of  the 


Monotone  Response  Theore™  for  Nonlinear.  Nonun'^o-" 
RC  Line's .  Given  a  nonlinear  RC  line  as  oescrioeo  aDove. 
Suppose  that  (because  of  circuit  parameter  uncertainty, 
the  use  of  linearized  models  for  nonlinear  elements, 
replacing  the  exact  input  by  input  bounds,  etc.,!  we 
do  one  re  of  the  following: 

a,  »-  -timate  the  input  e(t), 

b)  u>  estimate  one  or  more  R’s, 

c)  underestimate  one  or  more  C's. 


The  resulting  circuit  model  will  then  necessarily  over¬ 
estimate  the  output  vj(t)  at  each  instant  t  during  up 
transitions  (i.e.,  during  transitions  where  e  ^0, 
e  >  0  throughout. ) 

A  similar  result  holds  for  "down"  transitions  and 
estimate  errors  of  the  opposite  sign.  Using  part  a)  of 
the  assumptions,  this  theorem  allows  us  to  computationally 
propagate  upper  and  lower  signal  bounds  through  the  network. 
Using  parts  b)  and  c),  it  allows  us  to  replace  a  nonlinear 
line  0y  two  linear  ones,  one  strictly  faster  and  one 
strictly  slower,  to  which  the  linear  network  bounds  (4,5) 
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in  turn  apply.  We  nave  not  yet  succeeded  in  finding  a 
generalization  of  tms  result  that  will  apply  to  non- 
1  inear  RC  trees . 


III.  An  Apprpacn  to  Waveform  Bounding 
for  MQS  Logic  Gates 

Two  results  reported  nere  apply  to  HCS  uevice  models 
of  tne  form 

’d  1  r  f(VGB'  VDB’  vSB!'  (7) 


wnere  0,G,S  and  S  refer  to  drain,  gate,  source  and  sub- 
strate,  respectively.  For  specificity  we  consider  only 
n-cnannel  devices  m  this  paper.  No  special  algebraic 
form  for  f  is  assumed,  only  that  f  is  continuously 
di  fferentiaDle  and  satisfies  the  natural  monotomcity 
conditions 


—  -  0,  >  0,  <  0  (8) 

"VGB  3VD6  SB  *' 

evervwnere.  Thus  a  wide  variety  of  device  models  are 
allowed,  with  the  exception  that  (7)  does  not  allow  for 
snort-cnannel  effects. 


Our  aooroacn  will  be  to  reduce  a  mul ti pi e-i nput 
logic  gate  Dy  steps  to  an  equivalent  bounding  inverter' 
and  then  to  "find  bounds  for  the  response  of  this  inverter. 


3.11  Reduction  of  Series-Raral lei  Transistor  Networx 

tc  Equivalent  Bounoinc  Transistor 

We  nave  developed  a  metnod  for  reducing  any  series- 
parallel  transistor  networx  to  a  single  "equiva ’ ert 
bounding  transistor.  '  .sing  tne  tecnmgue  recursively, 
one  can  replace  tne  builup  or  pulldown  networx  oi  a 
multiple-input  cate  by  a  single  transistor  and  have 
--gnrous  bounds  for  tne  error  produced  oy  this  simpli¬ 
fication. 


For  example,  a  parallel  connection  of  N  transistors, 
all  identical  except  for  widths,  lengtns  and  gate 
voltages,  satisfies 

1  '■  *j(v3B’  VDB  ’  ’’SB '  & 

r  "v..,  ’  V3B ’  VS3'* 

.1-1  j  '  oo.. 

wnere  v-9  is  tne  vector  of  gate  voltages.  We  have  proven 
that,  because  of  tne_  assumptions  i 8 )_j_  there  exist  wea,  ^ 
-eo  independent  of  v-g,  ana  v-g  ana  vGn  tr.at  depend  on  v^g, 
sucn  that  .2  can  peJreolaced  oy  the  simpler  pounds 
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— JC 


'SB' 
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'or  all  vwg  _  vcg,  describing  a  single  transistor  with  a 
range  o*  gate  vot tapes.  Tne  ^unction  f  is  tne  same 
tnrougnout  ;9  ar.c  :10).  Figure  6  illustrates  this 
process  for  ■,  >  2 . 


Figure  6.  Replacing  a  parallel  transistor  network  By  an 


“equivalent  bounding  transistor  .  The  cost  of 
this  simplification  is  that  the  exact  value  of 
i  for  the  network  on  the  left  is  replaced  bv  a 
range  of  values  in  the  simpler  mooel  correspond¬ 
ing  to  V£g  i  vGB  1  vGg. 

3.2)  Reducinc  a  Multiple-Input  Gate  to  an  “Equivalent 

iSounninn  inverter1  ” 

A  gate  can  be  modelled  as  an  “equivalent  bounding 
inverter"  by  performing  the  reduction  outlined  in  section 
3.1  on  both  the  pull  up  and  pulldown  networks,  reducinc 
each  to  a  single  transistor.  Initial  trials,  comparing 
SPICE2  simulations  of  the  original  network  witn  simulations 
of  the  “eouivalent  boundino  inverter"  indicate  that  tne 
resulting  bounds  for  i0uI  ;vout)  differ  'rom  tie  exact 
valuet  by  only  about  +  TuS  for  orartical  circuits. 

3.3)  Bounding  the  Response  of  an  Inverte1-  and  Load  tc 
input  TransmonT 

When  applied  to  some  multiple-input  gates,  the  re¬ 
duction  procedure  described  in  the  previous  two  sub¬ 
sections  may  yield  an  inverter  in  whicn  the  oulluo  gate 
is  externally  driven.  But  for  simplicity  we  consider 
here  only  tne  case  of  a  standard  NMOS  depletion  -  load 
inverter  as  in  Fig.  7. 


figure  7.  Depletion-Ioaa  'nverter. 

To  bound  the  response  time  of  the  icaaea  inve^te' 
we  need  simple  bounds  on  the  function  i0ut  ( vpu- '  v'n 
wniro  is  the  difference  of  the  Dulluc  anaLpulicowr 
currents : 


I  X 


T  V  )  - 

out'  out’  in'  'pu''out'  'pd‘’out*  'in' 
Simple  linear  bounds  on  both  tne  Dulluc  anq  pulidowr 
currents  are  shown  in  Fie.  S.  The  resulting  ccurcs 
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Figure  8.  Simple  linear  bounds  on  the  duIIud  and  pull¬ 
down  currents.  The  latter  deoeno  on  V;  ,  and 
hence  on  t. 

Initial  simulations  using  this  approach  indicate  that 
the  delay  bounds  for  these  simplified  models  differ  f ra¬ 
the  delays  obtained  from  SPICE  simulations  Dy  about 
+  155. 

IV.  Further  Work  in  Progress 

Much  work  remains  to  be  done  before  the  theoretical 
basis  for  the  waveform  bounding  approach  to  timing 
analysis  is  complete.  Among  the  larger  remaining 
problems  are: 

1.  extending  the  Penfield-Rubinstein  bounds  to  in¬ 
corporate  time-varyinc  source  resistances,  such  as  those 
modelling  the  pulldown  current  in  Fig.  8, 
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2.  finding  bounds  for  the  response  of  an  RC  tree 
containing  pass  transistors, 

3.  investigating  the  tolerance  in  the  Pounds 
obtained  so  far  ana  fincing  tighter  ones  where  necessary, 
and 

4.  incorporating  effects  of  the  Mil ler  capacitance 
into  bounas. 
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